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Abstract 



0^ . We formulate Euler-Poincare and Lagrange-Poincare equations for systems with bro- 

ken symmetry. We specialize the general theory to present explicit equations of motion 
Q ■ for nematic systems, ranging from single nematic molecules to biaxial liquid crystals. The 

. geometric construction applies to order parameter spaces consisting of either unsigned unit 

pH I vectors (directors) or symmetric matrices (alignment tensors). On the Hamiltonian side, 

we provide the corresponding Poisson brackets in both Lie-Poisson and Hamilton-Poincare 
formulations. The explicit form of the helicity invariant for uniaxial nematics is also pre- 
sented, together with a whole class of invariant quantities (Casimirs) for two dimensional 
I incompressible flows. 

in 

^ : 1 Introduction 
o . 

^ ; 1.1 A simple example: the heavy top 

Symmetry breaking phenomena are widely common in several physical contexts, from classical 
^ . mechanics to particle physics. The simplest example is probably the heavy top dynamics, that is 
the motion of a rigid body with a fixed point in a gravitational field. This motion takes place on the 
special Euclidean group SE{3), which carries a semidirect product structure SE{3) = S0{3) (s)M^ 
involving the special orthogonal group SO (3). Semidirect product Lie group structures of this 
kind may be understood as "signatures" of broken symmetries in physical systems. For the case 
of the heavy top, the physical description involves a basic configuration space which is 50(3), 
i.e. the same as in rigid body dynamics. Although one can simply write Hamilton's equations on 
the cotangent bundle T*S0{3), one realizes that the heavy top dynamics is very different from 
the simple rigid body case. Indeed, the presence of gravity in the system produces a Hamiltonian 
which is not 5'0(3)-invariant, contrarily to what happens for the rigid body. Thus, besides the 
body angular momentum, one also considers the direction of gravity (in the body reference frame) 
as a dynamical variable in S"^. Upon identifying with unit vectors in M'^, one is leaded to enlarge 
the Lie group S0{3) thereby considering the semidirect product S'0(3) (s)]R^. Such an approach 
has a relevant infinite-dimensional analogue which yields the theory of compressible fluid flows. 
The main references for geometric mechanics on semidirect product Lie groups are [32, 22]. One 
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of the targets of the present paper is to generahze these semidirect product structures in the 
context of reduction by symmetry, thereby extending to the case when the vector space V in the 
semidirect product is replaced by a manifold M, so that the semidirect-product group 

structure breaks into the Cartesian product G x M. This introduction reviews the high points 
of the Euler-Poincare reduction theory and summarizes the results obtained in the paper, after a 
brief discussion of the geometric methods that will be used. 

1.2 Euler-Poincare approach to heavy tops 

This section briefly reviews the concrete mathematical background for the dynamics of the heavy 
top. This example will be still considered in the remainder of the paper for direct comparison. 

Although it is defined on the configuration space SO (3), the heavy top dynamics exhibits 
rotational symmetry breaking through the appearance of a potential term, depending on the 
direction of gravity es G S*^. The heavy top Lagrangian Leg : TS0{3) — t- M is written as 



where i3 C is the reference configuration of the body, p{A) is the mass density, M is the total 
mass, g is the acceleration of gravity, and ^ is the unit vector along the segment of length i that 
connects the fixed point of the body with its center of mass at t = 0. Evidently, the S0{3)- 
symmetry is broken: Leg{x, x) 7^ Le^{x~^x)- Rather, the system is invariant only with respect to 
S0{2) (rotations around the vertical axis), since this is the isotropy subgroup of e^. At this point, 
one denotes eg = Tq G and defines L : TS0{3)xR^ ^ M by L(x, X, To) := ^ro(x, x)- where Tq 
is interpreted as a new variable. We now let SO (3) act on the extended space TS0{3) x by the 
left action (x, X, Fq) H- {i^x^ i^X: "^ro), with the result that the function L is now S'0(3)-invariant 
since 



Thus, by S'0(3)-invariance, L induces the reduced Lagrangian / : so(3) x — M given by 




(1.1) 



MgiijTo ■ il^xC = MgiTo ■ xC, for all ^ E SO (3). 



where Cl = x '^X (so that 'hat' denotes the usual isomorphism so{3) ~ M'^) and F = x ^To- The 
moment of inertia tensor I is defined by (cf. [30]) 




Again, the reduced Lagrangian produces the Euler-Poincare variational principle 




to 



thereby yielding the equations 
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where the variations are evaluated as 6fl = S{x~^x) = S + 17 x S (with S := x~^^x) ^ind 
5r = O X r. This construction is referred to as Euler-Poincare approach for semidirect products 
(cf. [22]). In this particular case, the configuration Lie group is the well known special Euclidean 
group SE{3) := 5*0(3) (S)M^. Such a formulation is only possible upon identifying points on the 
sphere (e.g. the direction of gravity) with unit vectors in and by noticing that the modulus 
of the latter is preserved by the motion. Although in many situations this process leads to no 
contradiction, this is not always the case in condensed matter physics and it is useful to consider 
parameters (called 'order parameters') belonging to a manifold M rather than to a vector space 
V. For example, the heavy top would require a Lagrangian L : TS0{3) x S*^ — ?■ M. The next 
section gives an overview of the general setting in the context of condensed matter physics. 



1.3 Broken symmetries in condensed matter: nematic molecules 

An enlightening example of broken symmetry appears in the theory of nematic liquid crystals 
[6, 8]. In the simplest case of uniaxial molecules, such systems are continuum systems composed 
of oriented particles, i.e. particles that are endowed with a special orientation identified with an 
unsigned unit vector n ~ — n, called director. The presence of such a special direction in the 
system plays exactly the same role as in heavy top dynamics, as it is shown by the equations of 
motion for the single nematic particle: 

d 61 61 51 . 

— — = u X - — hnx— , n = vxn. (1.2) 

dtdu 6u on 

A direct comparison with heavy top dynamics shows the strict similarity between nematic particles 
and heavy tops. However, a director carries no sign and thus it does not belong to S^, but rather 
it takes values in 5'^/Z2, i.e. the unit sphere S"^ with diametrically opposite points identified, also 
known as the real projective plane MP^. Therefore the occurrence of the director breaks the full 
5*0(3) rotational symmetry of the particle, whose rotational motion is now invariant only under 
a subgroup V C 5*0(3), i.e. the isotropy subgroup of the director itself. For nematic particles, 
the isotropy subgroup of unsigned unit vectors is P = 0(2). (More precisely V = D^o, where the 
group Dao consists of rotations about the molecular axis and 180°-rotations about a normal to 
the molecular axis). In more general condensed matter applications, if O is the broken symmetry 
(the order parameter group) acting transitively on the order parameter space M, then the latter 
is isomorphic to the coset space O/V, where V is the isotropy subgroup V = Ono of a reference 
point no € M. In the case of nematics, the order parameter space is O/V = 50(3)/0(2) ^ MP^. 
Thus, after starting with the configuration space O (5*0(3) for nematics), the symmetry breaking 
requires to consider also dynamics on the order parameter space M ~ O/V. Notice that, this 
description refers only to the rotational motion of complex fluid particles and we have forgotten 
about its motion in physical space. However, the latter can be included a posteriori by considering 
Cartesian product configuration manifolds such a.s Q x O, where Q denotes physical space and 
we assume that O does not act on Q. In the first part of this paper we shall neglect motion in 
physical space in order to keep the treatment sufficiently compact. 

Nematic particles differ from the heavy top in that the order parameter space can be a man- 
ifold, rather than a vector space. Thus, the ordinary theory of heavy top dynamics needs to 
be extended for condensed matter applications. In particular, even if the order parameter coset 
space O/V is naturally associated to symmetry breaking, it is not clear how this coset space 
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arises from a fundamental approach in terms of reduction by symmetry. For example, the emer- 
gence of a broken symmetry often leads to consider O/V (rather than O) as the configuration 
manifold. In the context of nematic particles, the Ericksen-Leslie equations (cf. [27]) are simply 
the Euler-Lagrange equations on the projective plane RP^ ^ 5'0(3)/0(2). One of the tasks of 
this work is to show how this new configuration space emerges naturally from reduction theory. 
In more generality, 

this paper aims to provide a unified systematic framework for systems with broken symmetry by 
presenting the corresponding Euler-Poincare and Lagrange-Poincare approaches to formulate their 
dynamics. 

We shall mainly concentrate on the symmetry properties of finite dimensional systems, while 
the last part of the paper extends the treatment to consider the infinite dimensional cases. In 
terms of complex fluid systems, this means that we shall focus on the symmetry properties of 
the fluid single-particle dynamics. The question of how the underlying geometric structure is 
preserved in passing from the microscopic single-particle approach to the macroscopic continuum 
description will be the subject of future work. 

Notation. We have used the notation O for the order parameter group of a system with broken 
symmetry and have denoted by M the order parameter space, on which O act transitively. From 
now on we will use the notation O only when the group acts transitively on M . For general 
actions, the Lie group will be denoted by G. 

1.4 Geometric setting for symmetry breaking 

This section introduces some of the mathematical theory of systems with broken symmetry and 
anticipates how the coset space O/V emerges in a more rigorous framework. 

Let (9 be a Lie group (the order parameter group), acting transitively on the left of a manifold 
M (the order parameter space). Choose an element uq G M, and consider the isotropy subgroup 
V := OriQ- We have the isomorphism O jV ^ M, [x] = x'T^ ^ X^o? where V acts on O by right 
multiplication. The dynamics is described by a P-invariant Lagrangian L„q : TO — > M, which 
produces the Euler-Lagrange equations on TO. Now we notice that, from the transitivity of the 
action, any P-invariant function L^^ determines a unique function 

L -.TO ^ M L{v^, ipno) := Ln^iv^if) , with ip e O 

that is O-invariant under the right action oi il) E O given by {v-^.,n) )■ {v^il),i/j~^n). When 
is hyperregular, one can obtain the Hamiltonian formulation associated to a "P-invariant 
Hamiltonian '■ T*0 M by Legendre transformation. By proceeding analogously, one 
defines an (9-invariant Hamiltonian H{ay^,(pno) = Hngiaxf) '^^ T*0 x M (with the obvious 
notation eT*0). 

As explained in the previous section, this paper uses symmetry reduction theory to present 
the equations of motion on the corresponding reduced space. In particular, given a arbitrary left 
action n ^ gn oi a, Lie group G on a manifold M and a G-invariant function L : TG x M — > M 
under the right action {vh,n) i— )■ {vhg, g~^n), 

we develop two approaches: 
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(1) Euler-Poincare approach to symmetry breaking (Section 2). We make use of the 
diffeomorphism 

{TG X M)/G ^ X M, [vg, n] ^ {vgg~\ gn) , 

so that the function L induces a reduced Lagrangian / : g x M — j- M, which is defined by 
l{vgg~^ , gn) = L{vg,n). We obtain the equation of motion on the reduced space g x M. This 
construction generahzes the well known Euler-Poincare theory for semidirect products. We ex- 
plore the Hamiltonian side and determine the associated noncanonical Poisson brackets. In this 
context, we obtain a restriction of the Lie- Poisson bracket of [26], who considered the case when 
the manifold M is Poisson. Indeed, the Legendre transform of the Euler-Poincare equations yields 
the bracket in [2G] in the special case when M is endowed with the trivial Poisson structure. As 
from the discussion above, a particularly interesting situation is the case of a transitive action of 
G = O on M, so that M is isomorphic to the coset space O/V, V being the isotropy subgroup 
of O for a fixed element no G M . This result justifies the usual emergence of coset spaces for 
symmetry breaking in terms of Euler-Poincare reduction. 

(2) Lagrange-Poincare approach to symmetry breaking (Section 4). We shall derive 
Lagrange-Poincare equations by applying standard Lagrangian reduction to the ordinary La- 
grangian LriQ ■ TG — !■ M, defined by Lnoivg) := L{vg,no) for a fixed reference element no G M. 
Since L^^ is invariant under the isotropy group Gng of no, this process involves the quotient 
TG/Gno- simplicity, we suppose here that G = O acts transitively on M and denote by 
V the isotropy group of no (the case of a general action is treated in Section 4.2). In order 
to explicitly write the reduced equations, one needs a connection on the right principal bundle 
O ^ O/V ^ M. Using this connection, the reduced tangent bundle TO/V can be identified 
with the vector bundle TM ®m P over M ^ O/V, where © denotes Whitney sum, p denotes the 
adjoint bundle {O x p)/V and gothic fonts denote Lie algebras of corresponding Lie groups, as 
usual. This is a general construction [5], which is here applied to the coset bundle O O/V. 

1.5 Summary of main results 

After slightly extending the ordinary Euler-Poincare theory to account for order parameter 
manifolds, this construction is applied to recover the equations of motion for nematic par- 
ticles, thereafter formulating new Euler-Poincare equations for biaxial nematic particles and 
l^-shaped molecules. As a further step, the Euler-Poincare equations for ordinary nematic 
and biaxial particles are also expressed in terms of the corresponding alignment tensors, e.g. 
Q = 1/2 (nn^ — 1/31) for nematics. 

The second important result concerns the application of the Lagrange-Poincare reduction 
method to the coset bundle G G/Gno- Indeed, this method is shown to be extremely powerful 
when combined with the definition of a mechanical connection on G — ?■ G/Gq. The result of 
this combination is the formulation of the nematic particle motion on the configuration space 
M = SO{3)/0{2), that allows to relate systematically the variables (n, h) G TM with the 
corresponding Euler-Poincare variables [v, n) G so(3) x M. This step requires particular care 
because of the appearance of the extra constant r = n ■ u, belonging to the commutative Lie 
algebra o(2). For example, setting r = neglects rotations about n, consistently with the rod- 
like nature of nematic particles, while this is not sensible for heavy top dynamics, which turns 
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out to be also explained by exactly the same procedure (upon replacing 0(2) with 50(2)). It is 
worth mentioning that the possibility of a non-zero quantity r = n ■ i/ for the heavy top does not 
represent a limit of the theory: rather it represents a strong point of this approach that enables 
to split Euler- Lagrange (EL) dynamics on T{G/Gno) — '^^'^ from Euler-Poincare dynamics on 
the Lie algebra Qno = so (2) (i.e. trivial motion, in this specific case). As a final result of this 
method, we obtain the Ericksen-Leslie equation as a covariant EL equation on TS"^ and we extend 
its validity to the heavy top case, including a non-zero r G so (2). 

In the last part of this paper, the finite-dimensional treatment is extended to the infinite- 
dimensional fluid theory of liquid crystals. This is performed in two reduction stages: the first 
corresponding to the broken symmetry in the micromotion of nematic particles and the second 
corresponding to the fluid relabeling symmetry. On the Lagrange-Poincare side, this requires the 
use of a recent reduction procedure (known as 'metamorphosis'), first appeared in the study of 
the shape evolution in image dynamics [21]. After establishing the direct correspondence between 
the dynamics obtained by this approach and the Euler-Poincare reduced system, we perform 
the Legendre transform and apply a well known theorem [26] to produce two equivalent Poisson 
brackets for the same set of equations. Moreover, we present the explicit form of the helicity 
invariant for nematic liquid crystals, as well as a whole class of Casimirs functionals for two 
dimensional incompressible flows. 

2 Euler-Poincare equations and Lie-Poisson brackets 
2.1 The Euler-Poincare approach 

Let G be a Lie group acting on the left on a manifold M. We denote hj n ^-^ gn the action of 
g E G on n E M. Then G acts naturally on the right on TG x M via the free action 

{vg,n) ^ (vgh, h'^n) , 

where the action on the first factor is given by tangent lift of right translation on G, that is we 
denote Vgh := TRh{vg), where TRh is the tangent map to the right translation by the element 
h E G. Similarly, we will denote by agh the cotangent lifted action of right translation by an 
element h E G, that is, we have agh := T*Rh-i{ag) where evidently ag E T*G. The quotient 
space {TG x M) / G can be identified with g x M via the diffeomorphism 

[vg,n] ^ {vgg~^,gn). (2.1) 

• Assume that we have a function L : TG x Af — )■ M which is right G-invariant. 

• In particular, if uq E M, define the Lagrangian 

Ln, -.TG^R, Ln,{vg) := L{vg,no). (2.2) 

Then L„q is right invariant under the lift to TG of the right action of Gno on G, where Gng 
is the isotropy group of no. 

• Right G-invariance of L allows us to define / : g x M — )■ M by 

Kvg9~^, gn) = L{vg, n). (2.3) 
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• For a curve gt G G, let '■= QtQt ^ define the curve G M as the unique solution of 
the following differential equation with time dependent coefficients 

rit = (6)a/K), 

with initial condition no- Here (^t)Af G X(M) denotes the infinitesimal generator associated 
to the time dependent Lie algebra element C,t G g. The solution of this differential equation 
can be written as rit = QiUq. 

Theorem 2.1 With the preceding notations, the following are equivalent: 
(i) With Uq E M held fixed, Hamilton's variational principle 

ft2 



6 [\no{9,9)dt = 0, (2.4) 
Jti 



holds, for variations 6g of g vanishing at the endpoints. 
(ii) g satisfies the Euler- Lagrange equations for L„q on G. 
(iii) The constrained variational principle 



r"i2 

6 I l{^,n)dt = 0, (2.5) 

tl 



holds on g X M, upon using variations of the form 

where rj E g vanishes at the endpoints. 
(iv) The Euler- Poincare equations hold on g x M : 



61 61 (61 



dt6i ^ 6i \6n 

where J : T*M — g* , given by (J(am),.^) = (am, CM('Ti)), is the momentum map associated 
to the cotangent lifted action of G on T*M. 

The proof is an immediate generalization to that given in [22], in the case where M is a vector 
space on which G acts by a representation. 



Remark 2.2 The Euler-Poincare equations (2.6) can be rewritten in the form 

61 \ 
6^) 



d_ 

dt 



Ad*^ 



J ( 9~'f] ■ (2.7) 



Note that we have g G T*^M, since uq = g ^n. 
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Remark 2.3 (Relation to previous Euler-Poincare formulations) In the special case when 
M is a vector space acted on by ordinary representations, we recover the formulation in [22], use- 
ful for heavy tops and compressible fluids. Also, when the vector space M is acted on by afiine 
representations, the above picture recovers the results in [13], useful for complex fluids. 

Remark 2.4 (Terminology) Throughout the paper, we will refer to equations (2.6) as Euler- 
Poincare equations or Euler-Poincare equations for symmetry breaking since they generalize the 
[pure) Euler-Poincare equations 

d 61 61 ^ 

on a Lie algebra (see [30]). 

Remark 2.5 Since the Lagrangian L„q : TG — M is GnQ-invariant, it induces a reduced La- 
grangian tno '■ TG/G 

no ~^ Denoting by z„p : TG — ?■ TG x M the injection Vg i— )■ {vg,no), we 
can write L„q = L o By passing to the quotient spaces, the injection ino induces an injection 



: TG/Gno ^ X M, [vg]G^^^ ^ {vgQ \gno) 



whose image is g x Orb(?7,o), where Orb(?7,o) C M denotes the orbit of uq. We have the relation 
= I o between the two reduced Lagrangians ino ^■ 

Remark 2.6 (Including translational motion in physical space) As we remarked in the 
introduction, the present treatment neglects translational motion in physical space Q. It is 
important to emphasize that this simplification is completely irrelevant. Indeed, translational 
motion may always be included a posteriori, upon extending the Lagrangian L : TG x M — )■ M to 
be defined on the Cartesian product TQ x {TG x M). Thus, by assuming that G acts trivially 
on (5, the new Lagrangian Lng{q,q, g, g) will produce the extra Euler-Lagrange equation on TQ, 
which is then accompanied by the Euler-Poincare equation on g x M. In conclusion, the reduced 
Lagrangian / : TQ x (g x M) produces the equations 

d_6l_6I__ d_6l_ .*6I__ f6l_ 

dlTq Tq~ ' dtTi'^^^Ti' \6^ 

where one has to remember the relation h = ^m('^)- 



2.2 The Lie-Poisson equations 

Consider a function H : T*G x M — )■ M right invariant under the G-action on T*G x M given by 
(ttg, n) I—)- {agh, h~^n). In particular, the function Hno '■ T*G — j- M defined by Hnoio^g) = H{ag, Uq) 
is invariant under the induced action of the isotropy subgroup Gno '■= {g ^ G \ Q^q = ^o} for any 
no € M. The reduced Hamiltonian h : g* x M — )■ M is defined by h{agg~^, gn) := H{ag,n). 

Theorem 2.7 Fix an element uq G M. For a G T*G and fi := ag~^ G g*, the following are 
equivalent: 

i The curve a satisfies Hamilton's equations for Hno T*G. 
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ii The curve (/i, n) ^ q* x M is a solution of the Lie-Poisson equations 



fi = — ad*6h fi — 3 



6n 



n 



6h\ 

— (n), n(0) = no- 



These equations are Hamiltonian relative to the Poisson bracket 

'5f\ /^f5g\ 6f 



SI Sg_ 



Sn J Sfi 



6n J ' 6fi 



on Q* X M. 

As on the Lagrangian side, the evolution of the variable n is given by n = guQ. 



(2. 



(2.9) 



Proof. Canonical Hamilton's equations for Hng on T*G are equivalent to Hamilton's equation 
for H on T*G x M, endowed with the direct sum of the canonical Poisson bracket on T*G and the 
zero Poisson bracket on M, with initial value no- Since H is G invariant, one can apply Poisson 
reduction to obtain the reduced Hamiltonian equations on the quotient manifold {T*G x M)/G 
under the right action {ag,n) H- {agh, h~^n). We identify the quotient manifold (T*G x M)/G 
with Q* X M, via the diffeomorphism [agiTi] H- {agg~^ , gn). Using Proposition 2.1 in [26], in the 
particular case where the Poisson structure is trivial on M, we obtain the reduced bracket (2.9). 
One then observes that the Hamilton's equations associated to this bracket are given by (2.8). 



Legendre transformation. The preceding theorem is compatible with Theorem 2.1. Indeed, 
we can start with a Lagrangian L„q : TG — M as in §2.1, that is, we have a function L : TGxM — )■ 
M which is right G-invariant under the action [vg, n) i— )■ {vgh, h~^n), such that Lna{vg) = L{vg, no). 
Then L„,j is right invariant under the lift to TG of the right action of the isotropy group Gng 
on G. Suppose that the Legendre transformation FL„q is invertible and form the corresponding 
Hamiltonian Hno = Eno o ^L~^, where Eno is the energy of see [30]. Then the function 
H : T*G X M — !■ R so defined is G-invariant and one can apply this theorem. At the level of the 
reduced space, to a reduced Lagrangian / : g x M — )■ M we associate the reduced Hamiltonian 
/i : 0* X M — )■ R given by 

SI 

/i(/i,n) := (/i,^) - /(^,n), /i = 77. 

S^ 

Since 

Sh ^ ^ Sh SI 
— = C and — = - — , 
o/i on on 

we see that the Hamilton's equations for /i on g x M are equivalent to the Euler-Poincare equation 
(2.6) for / together with the kinematic equation 

n = ^M{n). 

Again, when M is a vector space acted on by a linear or affine representation, the above Lie- 
Poisson reduction recovers the results in [32] and [13], respectively. 
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3 Transitive actions: the dynamics of nematic molecules 



This section considers the interesting case of a Lie group G acting transitively from the left on 
a manifold M, called order parameter space. In this particular case we use the notation G = 
since we think of G as the order parameter group or the broken symmetry. 

This is a typical situation in condensed matter physics, which describes exactly the situation 
presented in the Introduction. The particular choice of the group O and the manifold M depends 
on the system under consideration. For example, the common case O = 5*0(3) and M = MP^ ^ 
5'^/Z2 yields ordinary nematic particles, while O = 5*0(3) and M = SO{3)/D2 yields biaxial 
nematics. Here D2 is the dihedral group generated by 180°-rotations and reflections. Many other 
choices are certainly possible and we address the reader to [33, 3 1] for reviews on these topics. 

Before proceeding, we emphasize that the present approach focuses on the rotational dynamics 
of a fixed nematic particle (e.g. ordinary or biaxial) in an external potential. Although this 
situation is physically improbable, its extension to more realistic configurations is straightforward. 
For example, translational motion in physical space Q can be included as discussed previously, 
while the dynamics of interacting particles requires extending the present treatment to the 
space x {O x O /V)^ . The continuum limit of a nematic lattice and the full hydrodynamic 
model are analyzed in section 5. 

As we have shown in the Introduction, if we fix a reference point rig € M and consider the 
isotropy group V = O„o = {ip E O \ tpriQ = rio}, then the orbit map 

C ^ M, X ^ Xno 

induces a diffeomorphism 

^no . ^ =xP^ xno . 

The reference point no corresponds to the coset [e] = eV = V oi O/V. Since the action is 
transitive, given a reference point no, a "P-invariant Lagrangian L^^ : TO — )■ M completely 
determines the right invariant function L : TO x [O/V) — ?■ M. Indeed, it suffices to define 
L{v^,n) := Lna{v^ip), where ip E O is such that n = ipno, by transitivity of the action. In 
particular, Ln^ determines the Lagrangians L^o for all other reference points mo G M. Note 
also that the injection ing is here a diffeormorphism since Orb(no) = M, that is, we have the 
diffeomorphism 

J^, : TO/V ^ X M, K] ^ {v^x~\ xno)- (3.1) 



3.1 Dynamics of uniaxial nematic molecules 

Theorem 2.1 applies to nematic particles without changes. Thus, in this case O = 50(3), O/V = 
RP^ and we have the following Euler-Poincare equations arising from the reduced Lagrangian 
I : 50(3) X MP2 ^ M: 

d 61 51 51 . ^, 

— — = X hnx— , n = iyxn (3.2) 

dtdu 5h> on 

where we have confused the space MP^ of unsigned unit vectors with ordinary unit vectors in M'^, 
thereby emphasizing the strict relation between rotational motion of nematic particles and the 
heavy top dynamics. Indeed, it is easier to consider the Lagrangian / : so(3) x [S'^jl^q) — M as 
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a function / : so(3) x S*^ — )• R possessing the Z2-symmetry /(^', v) = — v), for all v G 5o(3) 
and V G S*^. This is how equations (3.2) are written: one chooses a representative vector for n 
to compute the equations (i.e. a unit vector v G S'^ such that [v] = n), and then checks that the 
result does not depend on the chosen representative. 

The Lagrangian for the nematic particle is L^o '■ TS0{3) — )■ M 

1 2 

where j > is the moment of inertia and | ■ | denotes the norm associated to the S'0(3)-invariant 
Riemannian metric (X15X2) = Tr(x^X2) on S0{3). We notice that the inertia tensor I = jl 
leads to the schematic description of a nematic molecule in terms of a spherical rigid body 
endowed with a special direction uq. The unreduced potential energy is <l>no ■ SO (3) M. An 
example of such an external potential for a nematic particle is given by the quadratic expression: 
^no(x) = ^ |iio ■ /2, where A G M is a constant parameter and k plays the role of an 

external force field, e.g. an external magnetic (or electric) field. (Another possibility could 
be $no(x) = |no X (x~^k)|^/2). Besides the quadratic form of the potential, the important 
difference from the heavy-top case is that we now identify no with — ng and the same for the 
fixed vector k, which thus can point upwards or downwards without distinction. The presence of 
the potential $no breaks the symmetry, since Ln^ is only 0(2)-invariant (or better Doo-invariant, 
as mentioned in the Introduction) under right translation. If no G 5'^/Z2 is considered as an 
arbitrary director and we define L : TS0{3) x {S"^ /'L2) — )■ M by L(x, x,no) := Ln(,(x,x), then 
L is 5'0(3)-invariant under the right action (x, x, no) 1— )■ (x^, xV^, ^~^no). The argument is the 
same in the heavy-top CcLSG, clS shown by the following calculation: 

*(X, no) = ^ |no ■ X"^k|^ = ^ |xno ■ k|^ = $(xx~\ Xno) =: 0(xno) 

where we have introduced the reduced potential : 6*^/^2 — ?■ M. Thus, by S'0(3)-invariance, L 
induces the reduced Lagrangian / governing the rotational dynamics of a single nematic particle: 

K^,'^) = ^ J I^^P - </>(") = - ^ |n- k|^ (3.3) 

where v G so(3), and n G S'^/Z2. In this case, the Euler-Poincare equations for nematics reads 

jz> = V0(n) xn = A|n-k|kxn, n = i/xn (3.4) 

by (3.2). Moreover, the Legendre transform ^ = 6l/6u produces the reduced quadratic Hamilto- 
nian on so (3)* x MP^. 

^(M,n) = — + - |n-k|^ (3.5) 

The corresponding Lie-Poisson equations are 

d 6h 6h . 5h 

—H = — X fi+ — xn, n = — xn 

at ofj, on ofj, 

while the Lie-Poisson bracket of functions f,g^ ^-"(50(3)* x MP^) is 

{/, ,}(^, n) = X - J + n. ^- X - - - X - J . (3.6) 
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A special case of physical interest in microscopic theories of liquid crystals is the case of an 
ensemble of interacting nematic particles moving in the physical space Q. In this case, the 
unreduced Lagrangian L (o) (o) : TQ^ x TSO{?>)^ — )■ M is given by 

1 ^ 1 ^ 

i=l 1=1 iy^k 

where the first norm is associated to a Riemannian metric on Q, while O = 5*0(3)^ is endowed 
with the direct product group structure. The order parameter space is now M = {S'^/Z,2)^ ■ As 
before, Euler-Poincare reduction yields the reduced Lagrangian / : TQ^ x (so (3) x S'^/Z2)^ — )■ M, 
which is written as 

1 ^ 1 ^ 

/(g„i/„ni) = ^^Uif + ^JYI -^^{(li,(lk,T^i,nk) (3.7) 

i=l 1=1 i^k 

and produces Euler-Lagrange equations on TQ^ and Euler-Poincare equations on (so (3) x MP^)^. 
This approach is adopted in the physics literature to formulate microscopic approaches for nematic 
liquid crystals [35, 28, 2]. In this framework, when Q = M", the reduced potential is usually 
expressed as (cf. [28]) 

(j){qi, Qk, Ui, rifc) = V{qi - qu) + W{qi - qu) In^ -11^1^ 

where the angular factor |nj ■ n^l is reminiscent of the well known Maier-Saupe potential, in 
alternative to Onsager's expression |nj x n^l . In what follows, however, we shall pursue the 
question of how rotational motion is affected by symmetry breaking, without considering trans- 
lational motion in physical space, since the latter can be always considered a posteriori, once the 
reduction has been carried out on the rotational part of the Lagrangian. 

Remark 3.1 (The heavy top revisited) The same Euler-Poincare approach also applies to 
the simpler case when the order parameter manifold is M = S*^. In the case of symmetry 
breaking, the Lagrangian L : TS'0(3) — M is not S'0(3)-invariant. This is exactly the heavy-top 
case (see equation (1.1)), which involves a Lagrangian L that is only S'0(2)-invariant (rotations 
with respect to the vertical axis). The transitivity of the S'0(2)-action on 5*0(3) recovers the 
well known isomorphism S"^ = S0{3)/ S0{2), thereby producing a reduced Lagrangian / : so (3) x 
{SO{3)/SO{2)) —7- M. This approach recovers the well known treatment of heavy top dynamics 
(cf. e.g. [22]), avoiding the necessity of identify points in S"^ with vectors in M.^. 

3.2 Alignment tensor dynamics for uniaxial nematics 

Although the order parameter of a nematic liquid crystal often considered as a director, such a 
quantity is also difficult to work with analytically. In fact, the Landau-de Gennes theory typically 
involves the components of a symmetric matrix Q, usually known as alignment tensor. In matrix 
notation, this tensor is given by (cf. e.g. [6]) 

Q4(nn-4l) (3.8) 
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where nn^ is also denoted as the dyadic form nn. In order to write the Poisson bracket (3.6) in 
terms of the Q, we first need to compute functional derivatives. This is done by simply imposing 
the equality 5/(cr, n) = 5/(cr, Q), which yields the relation 



SI 



n 



SI 



At this point, it suffices to express the following term in terms of the alignment tensor: 



Tr 



Q ^ 



§^1 =Tr(Q 



SI Sg_ 
6fl ' 6Q 



where we have introduced the antisymmetric hat matrix fiij = Sijkfik and we have used the fact 
that Q is symmetric. Thus, the bracket (3.6) becomes 



{/,^}(M,n) = /X- 



S_l^Sg_ 



K 


'SI 


Sg' 




'Sg_ 


Sf 


)) 


6fi 


6Q 




5jl 


6Q 





(3.9) 



which is the Lie-Poisson structure associated to the semidirect product 5*0(3) (S) Sym(3), relative 
to the representation Q i— )■ xQx~^- Consequently, we have shown that the map (/x, n) i— (^t, Q) G 
50* (3) X Sym(3) is Poisson relative to the Lie-Poisson brackets (3.6) and (3.9). This result should 
not come as surprise, since the map n i— Q is equivariant with respect to the SO (3) actions 
n I—;- xn and Q i— )■ xQx~^- 

The resulting Lie-Poisson equations on so* (3) x Sym(3) are 



d 



6h 



' 6h' 



Q 



— Q 
6fl' 



with the notation Ai = e^kAjk and the following Hamiltonian (up to an irrelevant constant) 
arising from its previous expressions (3.5): 



M^,Q) = l|Mp + ^k^Qk. 



(3.10) 



Notice that more general potential terms of the type k^Q k2 are also allowed, which correspond 
to (n ■ ki) (n ■ k2) in (3.5). However, in the Landau-de Gennes theory, the potential energy 0(Q) 
involves more complicated expansions of the form 

0(Q)= Y^asnmcrW 

s,n>l 



where the coefficients asn are physical constants, see [8] . 

At this point, it is important to notice that the Euler-Poincare and Lie-Poisson equations 
involving the alignment tensor Q can also be derived directly by reduction of the Lagrangian 
I/Qo(x, x) or the Hamiltonian Hq^ix^ ^) (with (x, \E^) G T* SO {?>)), for a given Qo G Sym(3). This 
process is independent of the definition of the director variable n. 
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3.3 Dynamics of biaxial nematic molecules 

The same treatment applies to the case of biaxial nematic particles. In this context, we can 
identify the order parameter space with the manifold M of ordered couples n = (ni, of 
mutually orthogonal, unsigned unit vectors rii, 112 (cf. e.g. [29]). This choice is consistent with 
the microstructure of a biaxial molecule, which can thus be envisioned as a particle carrying 
two orthogonal director variables determining a special rotational state. The group 5*0(3) acts 
transitively on the ordered couples in M by matrix multiplication on each director: (ni,n2) 1— ?■ 
(xni,xn2). In order to determine the nature of symmetry breaking for biaxial particles, we can 
fix a reference couple no G M: for example, one considers ((1, 0, 0), (0, 1, 0)). Then, the associated 
isotropy subgroup SO (3) no C 5*0(3) is readily seen to be the dihedral group 

D2 = {diag(l, 1, 1), diag(-l, -1, 1), diag(-l, 1, -1), diag(l, -1, -1)} = 50(3)„„ . 

Therefore we can write M = SO{3)/D2 and we can express the Euler-Poincare equations for a 
biaxial nematic particle as 

d 61 51 51 . . , 

— — = 1/ X — + rii X — , ni = n,xu (3.11) 

atou oh> orij 

where we assume summation over repeated indexes in the first equation. A straightforward 
calculation may verify that the orthogonality condition ni ■ n2 = is consistently preserved at all 
times. A simple form of Lagrangian for a single biaxial nematic particle is given by 



|2 



/(z^,ni,n2) = ^ j \h'\^ - i(^Ai |ni ■ ki|^ + A2 |n2 ■ \^2? + A3 |ni x n2 ■ k2| 

where ki, k2 and ka again play the role of an external force field. A simple Legendre transform 
of equations (3.11) yields the following Poisson bracket for biaxial particles: 

r. 1. ^ f^f Sg\ f 5f 5g 5g 5f 

{/,,}(^,n„n2) - X +n,. X ^ - ^ X ^ 

fSf 6g 5g 6f\ 
+ n2- — X— - — X— . (3.12) 
V^n2 Ofi dn2 dfij 

The next section shows how the above bracket can be written in terms of the so called alignment 
tensors, analogously to our previous discussion concerning ordinary nematic molecules. 

Remark 3.2 (V-shaped molecules) The orthogonality condition rii ■ n2 = can be weakened 
by considering more general configurations of V^-shaped molecules carrying two directors spanning 
a fixed angle ^9. In this case, the isotropy subgroup is not the dihedral group D2. Rather, it 
is formed of vr- rotations only, and thus it can be identified with the cyclic group Z2. In this 
case, the order parameter manifold is M = 50(3)/Z2 and the equations (3.11) are still valid, 
thereby describing the dynamics of a single V^-shaped molecule with directors rii and 112, such 
that 111 ■ n2 = cos ^ at all times. 
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3.4 Alignment tensor dynamics for biaxial nematics 



In analogy with ordinary nematic molecules, we now show how one can write the Lie-Poisson 
equations for biaxial molecules in terms of the alignment tensor which is usually given by a linear 
combination of the two symmetric matrices 



- nin; 1 

2 V ^ 3 



B 



nsiio - Ini X ns) liii x ns 



(3.13) 



Following [37], we shall begin by considering the above matrices as two separate variables, so that 
the equality 5/(cr, rii, = Sf{(T, A, B) yields the relations 



6n2 



n 



112 X (ni X 112} 



5/ 
6B 



n,^ + n,x(nixn,) - 



In analogy to the procedure that we have followed for the case of ordinary nematic particles, we 
now express the following quantity in terms of (A, B): 



5/ ^9 ^ Sf 6g 
ni ■ — X h n2 ■ — X 



5^1 5\\2 



^ ■ rii X I rii — 



SI 



n2 X (ni X n2 



6B 
SB 



SI 
5^1 



Tr A 



ni X n 



SI 



T 



112 X T^2 



s_i 

6B 



, \ , \T Sf 

[ni X na) x (rii x 112) — 



5B 



S_l Sg_ 
6fi ' 6A 



\ f rSg Sf 

J-Tr(^n2n2-- 



+ Tr( (m X n2)(ni x n2f 





'SI 


Sg' 


^ +Tr(^B 


'Sf Sg- 




6fl 


6A 




_5fi ' 5B 



where we have used the same formulas that emerged in the case of nematic particles. Therefore, 
we can write the Poisson bracket (3.12) for biaxial particles in terms of the tensor order parameters 
as 



{/,^}(^,A,B) = Ai-(^x^)+Tr(A 



S± S^ 
_6fi ' 5k 



'Sg_ S_l 
6fi' 6k 



H 


'Sf Sg- 




'Sg_ 


Sf 


_6fi ' 6B 




5jji 


5B_ 



. (3.14) 



At this point, it is easy to see that any linear combination of the form Q = aA + /3B leads precisely 
to the bracket (3.9). Indeed, this follows upon noticing that 



S±^ S_l S_l^ 61 
Sk 6B ^(5Q 
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and by replacing the above functional derivatives in the bracket (3.14). 

Remark 3.3 (Other order parameter spaces) The present Euler-Poincare approach applies 
to any sort of order parameter manifold. For example, in superfluid dynamics one is faced to more 
complicated coset structures such as {S0{3) x 5*0(3)) / {S0{2) x Z2), which is the order param- 
eter space for the i?-phase of superfluid Helium-3 [31]. Another interesting example is provided 
by '^He-y4, i.e. the A-phase of '^He, whose order parameter space is (5*0(3) x 50(3)) /5'0(3) = 
5*0(3). In the context of liquid crystals one of the most complicated examples is provided by 
smectics, whose order parameter space involves the special Euclidean group [ > I]. Other possibili- 
ties still include more general complex fluids, such as micromorphic or micropolar fluids [10]. See 
[13] for a geometric treatment similar to the present one. 

Remark 3.4 (A fully symmetric case: magnetic moment dynamics) All the above ex- 
amples are indicative of how often symmetry breaking appears in physics. However, one should 
not forget that fully symmetric Lie-Poisson systems also emerge in condensed matter applications. 
A famous example is provided by the magnetic moments of ferromagnetic media, whose evolution 
takes place on coadjoint orbits of the rotation group 50(3). In this sense, the dynamics of the 
electron magnetic moment = 7/^ x H (where H is the external magnetic field and 7 is a physical 
constant) is a Lie-Poisson system on so*(3) with Hamiltonian h{iJ,) = 7^* -H. Thus, the magnetic 
moment dynamics possesses the same geometric interpretation as the rigid body dynamics, with 
the only difference that the latter also determines a geodesic flow on 5*0(3). A posteriori, one 
selects a special coadjoint orbit (the unit sphere 5"^) to be consistent with the definition of spin 
as a unit vector. Indeed, magnetic moments possess precisely the same symmetry properties of 
the usual spin variable, whose dynamics governs the theory of spin glasses [21]. 

4 Lagrange- Poincare approach to symmetry breaking 

While the last two sections presented the Euler-Poincare and Lie-Poisson formulations for systems 
with symmetry breaking, the present section will show how the same systems allow for another 
geometric description, whose underlying general theory is known under the name of Lagrangian 
reduction [3, 5]. While the Euler-Poincare theory applies to systems possessing a Lie group as 
the configuration space, Lagrangian reduction may be used to approach any system on a tangent 
bundle possessing a continuous symmetry. As we shall see, the resulting Lagrange-Poincare 
equations produce a special variable possessing a purely geometric character, and whose physical 
nature is analogue to the color charge in Yang-Mills theories. In these theories, the motion of a 
colored particle in a Yang-Mills field is a trajectory on a principal bundle B determined by an 
invariant Hamiltonian on T*B, where the Poisson bracket is canonical. In the particular case when 
the Hamiltonian is quadratic, the resulting Wong's geodesic equations produce the well known 
coadjoint dynamics for the color charge. It is clear that for Abelian Yang-Mills theories, the 
color is constant, like in the case of electromagnetism (cf. [30]). We shall see how a very similar 
situation also appears for the dynamics of nematic particles, whose celebrated Ericksen-Leslie 
equations are equivalent to the Euler-Lagrange equations on TMP^. 
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4.1 General Lagrangian reduction: a brief review 



Consider a right action ^iCxQ— >Qofa Lie group G on a manifold Q. Let L : TQ — )■ M be a 
G-invariant Lagrangian under the cotangent-lifted action of G on TQ. Because of this invariance, 
we get a well defined reduced Lagrangian £ : (TQ)/G — )■ R satisfying = L{vq). Assuming 

the group action is free and proper, the quotient space {TQ)/G is intrinsically a vector bundle 
over T{Q/G) with a fiber modeled on the Lie algebra q. Using a connection A on the principal 
bundle n : Q ^ Q/G we have a vector bundle isomorphism 



over Q/G, where the associated bundle g := Q Xq Q, is defined as the quotient space of Q x g 
relative to the right action (q,^) H- (<l>g(g), Ad^-i ^) of G. The elements of g are denoted by 
V = [q,i]G- Using the isomorphism a^, we can consider £ as a function defined on T{Q/G) © Q, 
and we write i{x, x, v) to emphasize the dependence of £ on (x, x) G T{Q /G) and v e g. However 
one should keep in mind that x, x, and v cannot be considered as being independent variables 
unless T{Q/G) and g are trivial bundles. 

We now formulate the Lagrangian reduction theorem. 

Theorem 4.1 The following conditions are equivalent: 
i Hamilton's variational principle 



aj, ■ {TQ)/G 



T(g/G)©0, K] 




to 



holds, for variations 6q{t) vanishing at the endpoints. 



ii The curve q{t) satisfies the Euler- Lagrange equations for L on TQ. 



iii The reduced variational principle 



S 




i{x, X, v)dt = 



holds, for variations 6x © of the curve x{t) ®v{t), where has the form 



V = —rj - [v, ri] + B{6x, x), 



(4.1) 



with the boundary conditions 5x{tj) = and rjiti) = 0, for i = 0,1. 



iv The following vertical and horizontal Lagrange- Poincare equations, hold: 




— {x, X, v) - -T7r(a:, x, v) = ( —{x, x, v), hB{x) 



> 



(4.2) 
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We now comment on the various expressions appearing in parts ii and iii. In the expression 
(4.1), D/Dt denotes the covariant time derivative of the curve f](t) G g associated to the principal 
connection A, that is, for f](t) = [q{t),^(t)]G, we have 

§-^[Qit),m]G = [lit), m + [Aiqit),qit)),m] 

The bracket [v, f]] denotes the Lie bracket induced by g on each fiber of g. The two-form B G 
Vt^{Q/G, g) is the curvature on the base Q/G induced by the curvature form E = dA + [A, A] G 
^^{Q, q) of Notice that for the formulation of the Lagrange-Poincare equations, the introduc- 
tion of an arbitrary connection V on the manifold Q/G is needed. For simplicity a torsion free 
connection is chosen. The partial derivatives 



di 

eT*{Q/G) and —{x,x,v) 



r- 7,* 
e 0X 



are the usual fiber derivatives of I in the vector bundles T{Q/G) and g, and 

^[x,x,v)eT:{Q/G) 

is the partial covariant derivative of I relative to the given connection V on Q/G and to the 
principal connection A on Q. We refer to [5] for details and proofs regarding the Lagrange- 
Poincare equations. Of course, there is an analogue result on the Hamiltonian side. Given a 
G invariant Hamiltonian H on T*Q, Poisson reduction yields the so called Hamilton- Poincare 
equations on the vector bundle T*{Q/G) © g*, see [1]. Given a Hamiltonian h = h(x, 7r,/i) : 
T*{Q/G) ©g*, the Hamilton-Poincare equations are 



D 
Dt 

X = 
dt' 



dl-L 



dh 

dn 



(4.3) 



-TT 



Vh 

dx 



/i, imiB{x] 



where dh/d7v and dh/dn are fiber derivatives, Vh/dx is the partial derivative defined in terms of 
the affine connection V on Q/G, see equations (16)-(19) in [4]. 



4.2 Lagrangian reduction on the coset bundle G — > G/Gn^ 

In this Section we specialize the situation described in §4.1 to consider the case of a Lie group 
configuration manifold Q = G, which is acted on by its isotropy subgroup where Hq is a fixed 
point on order parameter manifold M. By the quotient TG/Gng — T{G/Gno) ® flno (o^' 
Hamiltonian side: T*G /Gn^ — T* {G / G no)®Qno*) i this description provides the Lagrange-Poincare 
(or Hamilton-Poincare) approach to symmetry breaking, thereby extending the new phase space 
T*{G/Gno) to include a charge-like variable taking values in the dual isotropy subalgebra q^^. 
This process provides a rigorous geometric framework that naturally explains the emergence 
of the order parameter space G/Gn^ as the effective configuration manifold, starting from an 
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original system defined on the broken symmetry group G. This is the situation always appearing 
in physical applications (particularly in condensed matter physics) and its geometric formulation 
apparently differs from the Euler-Poincare approach, which does not treat G/Gno as the effective 
configuration space. 

In what follows, we suppose that the isotropy subgroup Gno C G is a Lie group, with 
dim Gno > 1- Since the G^g-action is free and proper, we have the right principal bundle tt^^ : 
G — >■ G/Gno — C)rb(no), g ^ gGno — gno- Given a principal connection A we have the usual 
vector bundle isomorphism 

TG/Gno ^ TOrb(no) © flTo, N ^ {iv99~')M{n), [g,Aivg)]Gj (4.4) 

over Orb(r2o), where n = gno and Qno is the adjoint bundle of the isotropy subalgebra Qno- Recall 
from Remark 2.5 that we have the diffeomorphism 

ino ■■ TG/Gno X Orb(no) C x M, [vg]G„^-, ^ {vgg~\gno). (4.5) 

Therefore, by composing with (4.4) we get the vector bundle isomorphism 

g X Orb(no) ^ T Orb(no) © s^o, i^,n) ^ (a^n), [g,Ai^g)]Gj (4.6) 

over Orb(no), where g & G is such that gno = n. 

As a concrete example of a principal connection, we consider the mechanical connection as- 
sociated to a Ad-invariant inner product 7 on g. The associated Riemannian metric on G is 
bi-invariant and is denoted by ((-, ■)). Recall that we have (see [31]) 

Aivg) = Iig)-\Jivg)). 

where the locked inertia tensor I{g) : g„„ g*^ is given by 

{i{g)vX) ■= {{VG{g),CGig))) = {{gv,90) = 7lg„^(^,C), with ^,r/ e g„o. 

Therefore the explicit expression of the locked inertia tensor is I{g)ri = 7]^^ (^7;-); where 7]^^ 
denotes the restriction of the inner product 7 on g to Qno C g. Note that when ^ G g then 
^{9)~^ -)\sno) ~ ^"o(0) where P.„o : g — )■ Qno is the orthogonal projector associated to 7 and 
7(^5-)lg,io denotes the restriction to Qno of the linear form 7(^,-) G g*. Indeed, one has 

I(^?)Pno(0 = 7L„(Pno(0,-) = 7(e,-)k 

The map J : TG — j- Q^o i^ given by 

{3{Vg)X) ■■= {{VgXGig))) = {{V9,g0) = 7(^^^^'a,C)L ' all C e Qno- 

Therefore 3{vg) = 7(5'~^'^5; -)|g„ aiid 

Aivg) = I{g)-\3ivg)) = I{g)-'[^ig-\,_)\^J = ¥noig~\) 
In this case, the vector bundle isomorphism (4.4) reads 

TG/Gno ^ TOrb(no) ©gTo, N ^ {ivgg-')Min), [g.fno{g-\)]Gj (4.7) 
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and the diffeomorphism (4.6) is 

X OrbK) ^ rOrb(no) ©sTo, {^,n) ^ (6,^, P„„(Adg-i Oh J (4-8) 

In order to compute the curvature of the mechanical connection, we use the formula B = dA + 
[A, A]. We have 

dA{ug, vg) = d{AiY))X - d{A{X))Y - Ai[X, Y]), 

where X, F G X{G) are two vector fields extending Ug,Vg. Using the left-invariant vector fields 
X{g) = gC, and Y{g) = gt], where ^, ?7 G g, we have 

dA{ug,Vg) = -Ai[X,Y]ig)) = -Aig[^,v]) = -^noi[^,v]) 

thus, we get 

SK,^;,) = [P„„(0,Pno(^)]-Pno([e,^]), ^ = 9"\, V = 9~\- (4.9) 

Suppose that we have chosen a fixed reference point uq G M. Using Lagrangian reduction 
(Theorem 4.1) for the G^g-invariant Lagrangian L^^ : TG — ?■ M and the reduced Lagrangian '■ 
T Orb(no) ©flno? we obtain the Lagrange- Poincare equations (4.2) in terms of [n, n) E T Orb(no)'. 

Remark 4.2 If the Lagrangian L^^ : TG — M is hyperregular, one can obtain the Hamiltonian 
description for the corresponding Hamiltonian : T*G — M. The reduction process is analogue 
and one ends up with the Hamilton- Poincare equations on the vector bundle T* Orb(no) © Qno* i 
see (4.3) with reduced Hamiltonian h„g for the general theory. We will refer to as the Lagrange- 
Poincare (LP) Lagrangian and to h„o as the Hamilton- Poincare (HP) Hamiltonian. 

4.3 Two equivalent approaches for symmetry breaking 

At this stage, it is useful to discuss the relation between the Euler-Poincare and Lagrange-Poincare 
descriptions associated to a G-invariant Lagrangian L : TG x M — > M. The equivalence of the two 
approaches arises as follows. Fix a reference point no G M and consider the induced Lagrangians 
: TG M, / : X M M, and ino '■ TOrb(no) ©orb(no) 0^ ^ I^- Then the following are 
equivalent: 

• (7 G G is a solution of the Euler- Lagrange equations for L^^. 

• ^ •= 99^^ ^ s-iid n := gn^ G M are solution of the Euler-Poincare equations (2.6) for /. 

• n := gno G Orb(rio) and ^ := [g,A{g)]G„^ G Qno are solutions of the Lagrange-Poincare 
equations (4.2) for 

Remark 4.3 (The case of transitive actions) In the case where G = (9 is an order parame- 
ter group acting transitively on the order parameter space M, we have Orb(no) = M, thus the 
reduced Lagrangian is defined on the vector bundle TM ®m P, where p is the Lie algebra of V. 
In particular the diffeomorphisms (4.4), (4.5), (4.6) become 

TO/V ^TM (BP, [v^]v ^ {iv^x~')Min), [x,Aiv^)]r) 
tn, : TO/V ^ X M, [n^]r ^ {v^x'\ xno) 
oxM^TM®p, {u,n)^ (i^a/ (n), [x, A{ux)]v) ■ 



20 



4.4 Lagrange-Poincare formulation of uniaxial nematics 

We now treat the particular case of nematics. Here the order parameter space is the projective 
plane M = RP^. Recall that MP^ is the non-orient able two-dimensional manifold given by the 
quotient of the two-sphere by the antipodal relation. An element n G MP^ is an equivalence class 
n = [v], where v G S*^ is a unit vector. The broken symmetry is the group of rotations 5*0(3) 
acting on the directors on the left by n i— )■ xn := [x^]- The isotropy group of a fixed direction no 
is the infinite dihedral group V = D^o generated by the rotations around the axis no and rotation 
by TT around an orthogonal axis. When the 2;-direction [(0, 0, 1)] is chosen for the reference axis, 
then the isotropy group D^o consists of matrices of the form 

The formula above shows how the group is actually isomorphic to 0(2), which is the group 
usually appearing in the condensed matter literature. 

Remark 4.4 The group i^oo should not be confused with Doo,h which contains also —I^ (i.e. 
minus the identity matrix) and is not a subgroup of 5*0(3) but of 0(3). Note that is not 
Abelian. For example the vr-rotation diag(l, — 1, — 1) does not commute with rotations around 
the vertical axis. 

The Lie algebra D^o of -Doo is given by matrices of the form 

-a \ 
a , (T G M. 
/ 

Using the hat map, it can be identified with the subspace Span(0, 0, 1) of on which the adjoint 
action acts by matrix multiplication. By identifying Span(0,0, 1) with the real line, and using 
the notations of (4.10) we obtain the adjoint action 

AdvK a = det(A^) a 

(recall that det{M) = ±1). The identification 50(3)/P'oo ^ ^P"^ is given by 

[X] G S0{3)/D^ ^n= [(xi3,X23,X33)] e MP^ ^ S^Z^. 

Giving a reference point no G MP^ and a Lagrangian L^^ : T5*0(3) — t- M describing the nematic 
particle, one can easily obtain the Euler-Poincare formulation associated to the Lagrangian / : 
so (3) X MP^ — M. In order to obtain the Lagrange-Poincare equations for i^o '■ TS0{3) / M, 
one needs to use the adjoint bundle. As we will see, it will be more comfortable to work with 
the sphere 5*^ instead of the projective plane. The adjoint bundle of the right principal bundle 
5*0(3) — )■ MP^ is the quotient space 

^ = (50(3) X Ooo) ID^ = (50(3) X M) /D^ 

relative to the right action of G D^o on (x,''") G 50(3) x M given by 

{X^r) ^ {xM,det{M)r). 
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Note that an element in the fiber (foo)^^ of the vector bundle Ooo MP^ reads [x,f], where 
X G S0{3) is such that n = xng, that is, [x] = n. The fact that also acts on r is due to the 
fact that Doo is not Abelian. 

In order to simplify the approach and to work with more explicit formulas, we replace the 
order parameter space MP^ with the two sphere S^. The breaking symmetry group 5*0(2) is now 
Abelian. The reduced Lagrangian is now defined on so (3) x S"^ but one has to recall that it is 
invariant under a change of sign for the variable in S"^. In this case, the isotropy group is 5*0(2), 

and the adjoint bundle so(2) = 50(3) y<so{2) M is a trivial bundle, since 50(2) is Abelian and 
thus the adjoint action on its Lie algebra so (2) ~ M is trivial. 

We now describe the mechanical connection associated to the Ad-invariant inner product 
7(^,77) = |trace(^^?7) on so(3). Recall that the Lie algebra so(3) and so(2) are identified with 

and M(0, 0, 1), via the hat map. On the inner product 7 is the standard inner product and 
the projection P : so(3) ~ — )• so(2) ~ M is simply given by taking the third component. The 
mechanical connection is thus given by 

Aiv^) = nx-\) = ix~\h- (4.11) 
According to formula (4.9), the curvature is 

B{U^, V^) = - {U X k)3 = -Z/i K2 + Z/2 = 1^ = X'^V^- 

In order to compute the reduced curvature on 5^ we first note that the tangent map to the 
projection TTno : 50(3) — )■ 5^, x ^ X^o reads 

TTTno : T50(3) TS^, />x ^^ X xno = X n. 

Note that if h G T^S^ is given, then we have TTi^^in x ri^) = n, where the 'hat' map denotes 
the usual isomorphism ^ so (3). Therefore, the reduced curvature is 

i3n(n, rii) = i3 x ri^, n x riixj = ~ (x ""^(n x ri) x x"^(n x m))^ 

= - ((no X x"^n) X (no x x~^i^))^ = - (x^ri x X^rii)3 

= — (x^(ri X rii))g = — n-(ri x m). (4-12) 

In the last equality we used the formula {x^i^)3 = Xa^i = n • i/, valid since x ^ 50(3) is such 
that xno = n, where no = (0,0, 1). Up to a sign the curvature is given by the volume of the 
polytope generated by the vectors n, ri and rii. 

Using the mechanical connection, and the fact that the adjoint bundle is trivial, the diffeo- 
morphism (4.8) reads 

so(3) X 5^ ^ r5^ X M, (i/,n) h4 (n,i/ X n,n-ry) = (n,ri,r). (4.13) 

Indeed, the infinitesimal generator associated to v E so(3) ^ M'^ reads as t'52(n) = i/ x n G Tn5^ 
and the second component of (4.8) is 

P(Ad^-i u) = (x^'i^)3 = (X^i^)3 = Xai^i = n-J^, 

In terms of equivalence classes [v^ G T50(3)/50(2), the variables n and u are given by n = xno, 
and V = v^x~^ ■ 
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We now obtain the formula for the inverse of the vector bundle map (4.13). Denoting r := n u 
and h := 1/ X n we have n x (i/ x n) = {n n)u — n(n-i/) = ly — rn, and 



i/ = nx(i/xn)+rn = nxn + rn. (4-14) 

This proves that the inverse of the vector bundle map (4.13) is 

(n, h, r) G TS*^ x M i-)- (n x n + rn, n) = (t/, n) G so(3) x S*^. 

Therefore, the Lagrangian / in (3.3) and the LP Lagrangian (now simply denoted by i, recall 
that no = (0, 0, 1)), are related by 

n) = i{n, u X n,n u), i{n, h, r) = /(n x ri + rn, n). (4-15) 

We now compute the Lagrange-Poincare equations associated to the reduced Lagrangian i = 
i{n, h, r) : TS^ x M — M. Since the group 5*0(2) is Abelian, the Lie bracket is zero, the covariant 
derivative coincides with usual derivative. Using formula (4.12) for the reduced curvature of the 
mechanical connection, the Lagrange-Poincare equations (4.2) read 

ddi_^ V d£ a 16) 

dt dr ' dt dh dn dr 

where the curvature term is evidently given by n x n = —1^3^, and V/dt, V£/dn are associated 
to the Levi-Civita connection on S"^ induced by the inner product on M^. The explicit form of 
the Lagrangian can be written immediately upon using (4.15) and recalling equation (3.7). 
We obtain the LP Lagrangian 

^(n, n, r) = /(n x ri + rn, n) = - j |n x ri + rnp — 0(n) 
= |n X rip + ^jr^ - 0(n) 

= ^J|n|2 + ijr2-0(n). (4.17) 

where the potential term may be expressed again by the quadratic expression 0(n) = A |n ■ k|^, 
or alternatively by 0(n) = A |n x k|^, for a fixed number A and an unsigned unit vector k. Upon 
inserting the expression (4.17) in the Lagrange-Poincare equations (4.16) we get 

. ^ .V. V</. . 
r = U, j—n+— — =jrnxn. 
dt on 

Note that does not depend on ri, so V0/(9n coincides with the gradient of the map (p : S"^ ^M. 
with respect to the Riemannian metric induced by M'^ on 5^. In order to recover the usual form 
of the equations, we will interpret n as a curve in and as a map defined on M^. In this case 
we get the equation 

j ri — 2gn + V0(n) = j r ri x n with 2g := n • (j ri + V0(n)) 

and where V0 denotes the gradient of view as a map on M^. This equation evidently reduces 
to the celebrated Ericksen-Leslie equation for the case r = (cf. [6]). 
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Remark 4.5 (Physical nature of the variable r) The conserved variable r = Uq ■ Uq is evi- 
dently the projection of the angular velocity of the nematic molecule on its director. Thus, when 
this quantity is non-zero, it encodes the effect of rotations of the molecule about the director. 
Because of the particular form of the Lagrangian, this effect is taken into account only by the 
curvature term on the right hand side of the Euler-Lagrange equation. An analogous situation 
holds for the heavy top dynamics (see remark 4.6 below). However, for nematic molecules, such 
rotations about the director are irrelevant, due to the rod-like nature of nematics. Thus, the 
convention r = producing Ericksen-Leslie dynamics is the most natural in this case. On the 
other hand, this is not true for the dynamics of the heavy top, which is a rigid body (with fixed 
point) of arbitrary shape. 

For the Hamiltonian side, one needs the dual vector bundle map of (4.15) over S"^. We have 

(/X, n) G so(3)* X ^ (n, tt, w) = (n, ^ x n, /x-n) G T*S^ x M*, (4.18) 

whose inverse is given by 

[n, TT, wj G T*S^ X R* (n X TT + wn, n) = (/x, n) G so(3)* x S^. 

Therefore, the Hamiltonian h in (3.5) produces the corresponding HP Hamiltonian by h(n, tt, w) = 
h{n X TT + wn, n). The nematic particle Hamiltonian reads as 

h(n,7r,uO = — IttP + — + 0(n) (4.19) 
2j 2j 

which evidently differs from the Lie-Poisson Hamiltonian 

^(/^>n) = + (f){n) 

only by the constant factor w, without changing the content in physical information. In the 
case of (4.19), the quadratic nature of the Hamiltonian produces a Kaluza-Klein construction 
analogue to that underlying Wong's equations in Yang-Mills theory. In this setting, substitution 
of the Hamiltonian (4.19) in the the Hamilton-Poincare equations (4.3) yields 

. „ .. V V0 u; 
w = 0, jn = 77, — TT = — ttt X n, 

at on J 

so that the usual Ericksen-Leslie equation is recovered for w = 0. 

Remark 4.6 (Comparison with the heavy top) It is evident that the above reduction can 
be performed equivalently for the heavy top dynamics by expanding the quotient TS0{3) / 5*0(2) ~ 
TS^ © M, which gives analogous equations to those above. Although the Lagrangian (and the 
Hamiltonian) keep the same form as the single nematic particle, the main difference resides in 
the form of the potential, whose explicit expression is now (f){T) = A F ■ k, where we use the 
notation of §1.2. This expression is evidently 5'0(2)-invariant, rather than Doo-invariant, consis- 
tently with the broken symmetry of the heavy top system. Indeed, in the heavy top case, the 
dynamical variable n belongs to the sphere S*^ ~ S0{3)/ S0{2) rather than to the projective 
plane MP^ ~ S0{3)/Doo, thereby reflecting the different nature of the two breaking symmetry 
subgroups 50(2) and D^. Moreover, we remark that the Ericksen-Leslie equation also holds for 
the heavy top dynamics provided the potential is now (f){T) = A F ■ k. However, in this case, 
setting r = has no physical motivation. 
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Remark 4.7 (Lagrange-Poincare approach for moving nematic particles) From a phys- 
ical point of view, one should also take into account the translational motion of nematic particles 
in the physical space Q, thereby extending the LP Lagrangian to I : TQ x {TS'^ x M) — M. 
This process has the only effect of producing an extra Euler-Lagrange equation on physical space, 
without changing any of the present geometric construction. The final equations are then 

±dl_m_^ Z^_Z^-^nxri ——-0 
dt dq dq ' dt dh dn dr ' dt dr 

with the Lagrangian 

^(g,g,n, n,r) = ^ \\qf + ^ j \h\'^ - $(n) + ^ jr^ 

Remark 4.8 (Biaxial nematics) Recall that the case of biaxial nematic particles involves the 
order parameter space S0{3)/ D2, where D2 is the finite dihedral group. This discrete symme- 
try cannot be considered under the preceding Lagrange-Poincare approach, since the latter is 
defined only for order parameter spaces O/V, involving an isotropy subgroup V <Z O, such that 
dim(7') > 1. When the isotropy subgroup V is discrete, then the situation requires more care in 
order to take into account the trivial nature of its Lie algebra p = {0}. The broken symmetry 
group O becomes a principal bundle O ^ O/V with discrete fiber and the reduction process in 
this case remains unknown. 



4.5 Summary for uniaxial nematics 

This section gives an overview that summarizes the different approaches that have been carried 
out so far for uniaxial nematic particles. The starting point is the the unreduced S'0(3)-invariant 
Lagrangian L : TS0{3) x S"^ — ?■ M that describes the dynamics of a single molecule. Then, we fix 
a reference direction uq and consider the induced Lagrangians 

Ln„ : TS0{3) R, Ln„ = L„„ (v^) 

l:30{3)xS^ — >R, l = l{u,n) 
£„o -.TS^ xR-^R, = 4o (n, ri, r ) , 

where we choose the reference configuration no = (0, 0, 1) for simplicity. Given a curve x ^ S0{3) 
and the vectors n = x^o ^ S'^ and v = XX~^ ^ so (3), the following are equivalent: 

i The curve x is a solution of the Euler-Lagrange equations for Lng- 

ii The curves v and n are solutions of the Euler-Poincare equations for I: 

d 61 51 61 . ,^ , 

— - h — xi/ = nx— n = i/xn. (4.20) 

dt6u 6v 6n 



iii The curves n and r are solutions of the Lagrange-Poincare equations for 



-no- 



d dC d df df df 

dt dr ' dt dh dn dr ' ' ' 
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To obtain the equations (4.20) from the Euler-Poincare equations (2.6), we used the formulas 
ad* K = —V X K and the expression of the momentum map 

J : T*5^ ^so(3)*, J(n,7r) = n X TT. 

for the cotangent hfted action of S'0(3) on T*S^. The conservation law (2.7) for nematics reads 



d f ^^61 



1 61 

On can pass from the Euler-Poincare (4.20) to the Lagrange-Poincare equations (4.21), using the 
relations: 

1/ = n X n + rn 
r = nv, h = v X n. 

The link between the Lagrangian variables (x, x) ^-iid the reduced variables {u, n) and (n, ri, r) is 
illustrated in the following diagram. 



TS0{3) 

(x,x) 



TSO{3)/SO{2) 

[{x,x)] 




TS^ X R 

(xno, XX"^ X Xno, XX~^ ' X'^o) =■ (n, n, r) 



Remark 4.9 (unit sphere vs. unit vectors) Instead of the sphere M = S'^, one can start 
with the vector space M = M'^ on which 5*0(3) acts by matrix multiplication. This is the point 
of view adopted in [13]. In this case, if one chooses the initial condition no G S"^, then we still 
have n(t) = x(^)no G S"^, thus these approaches are equivalent. 

Remark 4.10 (Directors and the real projective plane) Recall that the parameter space 
for nematics is the projective plane MP^ and not the sphere S"^. However, the results i — iv are 
still true when one works with the projective plane, that is, when n is interpreted as a director 
instead as a unit vector. The only change concerns the Lagrange-Poincare equations which are 
less explicit in the case of the projective plane. This is the reason why we use S"^ instead of MP^ 
above. 

Remark 4.11 (Euler-Poincare equations via Lagrangian reduction) From our discussions 
we realize that the Euler-Poincare and Lagrange-Poincare equations arise from two different ap- 
proaches in reduction theory. However it is possible to obtain the Euler-Poincare equations by 
standard Lagrangian reduction, provided one suitably enlarges the physical configuration space. 
We address the reader to [3, 11] for the case of order parameter vector spaces. This picture 
produces the so called Clebsch constrained variational principle [JH]. 
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5 Hydrodynamics of nematic liquid crystals 



While the previous discussions have focused on the geometric dynamics of a single particle with 
broken symmetry, this section develops the same ideas in the more physical situation of a fluid 
system of particles with micro-structure, such as nematic particles. Thus, this section provides 
the link between the geometric treatment previously applied at the microscopic single-particle 
level and the macroscopic models usually adopted for liquid crystal dynamics. We shall see how 
the same ideas apply to continuum media without substantial modifications. The extension to 
nematic fluid dynamics requires considering the diffeomorphism group, i.e. the particle relabeling 
group well known in the Lagrangian picture of fluid dynamics. 

We begin by recalling our notation: O be the order parameter group of a certain particle with 
broken symmetry, and denote by M the order parameter space, on which O acts transitively, with 
isotropy group V. As we have seen before, passing from one particle to a system of particles 
simply consists in replacing the group O by the direct product O'^ acting on the cartesian product 
M^. When a continuum of particles with broken symmetry is considered, one needs to consider 
as symmetry group the group G := J^{T>, O) of smooth maps form the physical space D to the 
order parameter group O of the single particle. We still denote \yy x '■ ^ ^ these maps. The 
order parameter space is the manifold of maps M. := J^{T>, M) on which G acts by the naturally 
induced pointwise action. 

In order to describe the hydrodynamic of such systems of particles (called complex fluids), 
we need to include the group of diffeomorphisms in the symmetry group. For simplicity, we 
suppose that the fluid is incompressible, since the compressible case would require only a slight 
modification. Thus we need to consider the group Diffvoi(^^) of volume preserving diffeomorphisms 
of relative to a fixed volume form /i on V. It will be convenient to fix a Riemannian metric g 
on V and to choose /i as the volume form associated to the metric. 

5.1 Euler-Poincare formulation 

In the Euler-Poincare framework, the dynamics of complex fluids is obtained by considering the 
semidirect product G = DiS^oi(J^)®J^{1^,0) acting on advected variables by linear or affine 
representation as shown in [13]. Here the advected variable is the order parameter function 
n e J^{T>, M), on which G acts by the left action 



Recall that denotes the action of J^(T>, O) on J^{T>^ M) naturally induced by the action 
of O on M. One easily checks that (5.1) defines an action of the semidirect product G = 
Difr™i(r') (s) J^(P, O). Using the expression 



n ^ {r],x)n := (x^) or] 



-1 



ij].x) e G. 



(5.1) 



ad(„ j,)(m, k) = {£ufn + k- Vz/, £ui^ + ad* k) , 



for the infinitesimal coadjoint action, we obtain from (2.6) the equations 




div u = 



(5.2) 
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together with the kinematic equation 

n + u-Vn = vm ° 

The Lagrangian Lr^,,:T{ Diff™i(r') (s) O)) ^ M is usually of the form 

Jv ^ Jv JV 

where the norm of fj is given by the Riemannian metric g, the norm of x is associated to a right- 
invariant Riemannian metric on (9, and F is the free energy. One easily checks that this expression 
is invariant under the right action of (?7, x) G G = Diffvoi('C) (§)-7-'(P, (9) given by cotangent lift 
on TG and acting on uq as uq v-^ X~^('^o ° ^7)- The reduced Euler-Poincare Lagrangian recovers 
the expression 

l{u,v,n) = - I / / F{n,\/n)fi, (5.3) 



2 Jti 2 
where F is the free energy. 

5.1.1 Euler-Poincare fluid equations for nematic molecules 

The case of nematic liquid crystals can be treated as above, choosing O = S0{3), M = RP^. In 
this case the free energy F is usually given by the Oseen-Zocher-Frank expression 



F(n, Vn) = ^iri(divn)^ + ^i^2(n ■ curln)^ + ^-ft'sln x curing, 



(5.4) 



where the constants Ki, K2, K3 are respectively associated to the three principal distinct director 
axis deformations in nematics, namely, splay, twist, and bend. At this point, the system (5.2) 
produces the equations 



dfU + n- Vn = —di { Vn 



dF 



Vp, 


div u 


= 


SI 


9^ ^ 


fdF 


dn 




an 





(5.5) 



j iy + u-S/v) = h X n, h := 
together with the kinematic equation 

h + M- Vn = 1/ X n. 

The Ericksen- Leslie fluid equations follow easily from the above relations, under the assumption 
that the initial condition verifies i^o-uq = [13]. A direct computation using (2.9) shows that the 
associated Poisson bracket reads 



{/,^}("^,M,n) 



5/ ^ 
(5m' 5m 



^ (^x^ + V^-^-V^-^ 
(5/^ (5/x (5/^ 5m 5n 5m 



5f 5f 5g 

n X — + Vn - — ,— 
(5/x 5m 5n 



5g 5g 5f 

nx — + Vn-— , — 
5fx 5m on 



(5.6) 
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where m = 61/ 6u is the fluid momentum and the brackets ( , ) denote duahty. The first two 
terms is the Lie-Poisson bracket, associated to the semidirect product Diffvoi(I^) ®J^(T>, 5*0(3)). 
Of course, this bracket is consistent with that of the single particle, see (3.6). 

From the discussion above, it is easy to generalize the Euler-Poincare fluid equations for 
uniaxial nematics to the case of biaxial liquid crystals. Indeed, upon choosing O = 5*0(3) and 
M = 5*0(3)7-02, these equations can be written down directly as follows: 



d 61 
dt 6u 

d 61 
di~6^ 



61 


61 


6u 


6^ 


[61 













61 



5^Vn.---Vp 



1=1 



6l_ 
6^ 



1=1 



6ni 
6ni 



div u = 
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together with the kinematic equation 

ri; + n- Vn^ = X 1/ , / = 1, 2. 
On the Hamiltonian side, the corresponding Poisson bracket reads 



{/,fi'}("^,M,ni,n2) 



5/ 5^ 
5m' 6m 



6^1 6fi 6fi 6m 6fj, 6m 



1=1 



51 
6^1 



51 5g_ 
6m^ 6ni 



6g 6g 5f 

n; X — + Vni ■ — , — 
o/j, dm on/ 



The next section extends the Euler-Poincare approach for the alignment tensor dynamics of a 
single nematic molecule to the case of nematic liquid crystals. 



5.1.2 Alignment tensor dynamics in nematic liquid crystals 

It is well known that fluid equations are obtained from the single particle dynamics by taking 
statistical averages with respect to some probability distribution function. The case of liquid 
crystals is not an exception to this procedure and one is usually interested in the dynamics of 
the averaged alignment tensor Q, where the 'bar' symbol denotes an appropriate averaging of the 
single particle alignment tensor Q. In order to formulate the Euler-Poincare equations for liquid 
crystals in terms of the averaged quantity Q, one proceeds by simply replacing the order parameter 
space M = MP^ with the space of symmetric matrices Sym(3). The Euler-Poincare reduction 
process is performed on the unreduced Lagrangian Lq„ : T{DiS^o\(P)®J^{V, S0{3))) ^ M with 
respect to the action 

where evidently Q G J^(T>, Sym(3)) and we have suppressed the averaging notation. Thus, if one 
specializes the general equations (5.2) to the case under consideration, one obtains 



' d 61 ,61 61 ^ 

+ £u— + — -Vi^ = -VQ 



dt 6u 
d 61 

+ div I —u 
at 6u V 6u 



6u 6^ 

61 



61 



V X 



6l_ 
6^ 



Q, 



61_ 
6Q 



div n = 
Q + n-VQ = [i>,Q] 



(5.7) 
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or more explicitly, for a Lagrangian of the form (5.3) with n replaced by Q, 



dtu + u-Vu = —dk 



dF 



j{u + u-Vu) = [Q,H 
where the molecular field H is given by 



divM = 



Q + n-VQ = [i>,Q] 



61 dF 



^--6Q-dQ'^' 



dF 



with the notation Ai = EijkAjk and the summation convention over repeated tensor indexes. 
The free energy F{Q, VQ) appearing in the reduced Lagrangian l{u, i/, Q) is usually given by the 
Landau-de Gennes free energy in the form [7] 

F(Q, VQ) = Ki II VQf + K2\\V -Qf + K3 Tr(Q V x Q) + aai Tr(Q2) + 031 Tt{Q^) 



where the tensor norms are given by total contraction of the indexes, i.e. || VQ||^ = XI (^feQ 
and II V ■ Q||^ = ^ (^diQ^WkQ^-^) , while (Vx)^^. = eijkdk is considered as a matrix operator. Higher 
order expansions in the alignment tensor are also possible. 

On the Hamiltonian side, the corresponding Poisson bracket reads 



{/,^}("^,M,Q) 



+ Tr 



5/ 5^ 
(5m' 6m 



Q, 



6fi 



5i 
6Q 



— X — - 

6n 6fx 

6m J 6Q 



.5f 6g 



■ 5g 6f 



6n 6m 6fi 6m 



Tr 



Q, 



Si 

6fi 



6Q 



+ 



6m J 6Q 



and it is easy to recognize that the same relations also hold for biaxial liquid crystals. 

Remark 5.1 (Compressible fluid flows) In the Euler-Poincare setting, the generalization to 
compressible fluid flows is straightforward. Indeed, it is sufficient to consider the whole diffeo- 
morphism group Diff(X') and to enlarge the order parameter space by replacing J^{V,M) with 
the product J^(T>,M) x Den(D), where Den(D) denotes the space of densities on V. Then one 
repeats the construction above, upon considering the following action of G = Diff (P) (§)J^{V, O) 
on J^{V,M) X Den(P) 

where J(?7) is the Jacobian of rj with respect to a volume form on V. In this sense, compressibility 
is by itself another example of symmetry breaking in infinite dimensions (cf. [22]). Indeed, even 
in the case of ordinary isotropic compressible fluids, the full Diff-symmetry of the unreduced 
Lagrangian : TDiff(D) — )■ M is broken by the presence of the density variable po ^ Den(D), 
so that the only symmetry which is left is given by the infinite-dimensional isotropy subgroup 
DiffpQ(P). Upon applying the Euler-Poincare approach, one performs the reduction in (3.1): 
TDiff('D)/Diffp(,(r') ~ X{V) x Den(P) because of the transitivity of the action (see e.g. [25]). 
In principle, it is also possible to apply Lagrangian reduction and write the Lagrange-Poincare 

equations on the reduced space TDiff(r')/Diffp(,(r') ~ TDen(r') © X^oiiV). Although the next 
discussion is devoted to Lagrangian reduction, we shall consider again the case of incompressible 
nematic liquid crystals. 
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5.2 Lagrange-Poincare formulation of uniaxial nematic liquid crystals 



The Lagrange-Poincare approach to continuum nematic media produces another variant of the 
fluid equations for hquid crystals. One starts with the same Lagrangian as before, namely 

Lno{v,V,X,x) = 7: I + ^ / jIxIV- / ^ ((xno) or/-\ V((xno) or/"^)) /i. 

Jv Jv Jv 

However we will now reduce the dynamics in two steps. The first step is a Lagrange-Poincare 
reduction involving the internal variables. The second step involves the diffeomorphism group 
and is done by using the metamorphosis approach. 

5.2.1 First reduction stage: Lagrange-Poincare approach 

As a first step of reduction, one applies the Lagrange-Poincare approach for Ln,, according to the 
vector bundle isomorphism 

TBiS^oiiV) xTJ^{V,SO{3))/J^{V,SO{2)) TBiS,oi{V) x {TT{V, S^) ® T{V)) 

over Diffvoi('C') x J^(V, S"^). By a slight abuse of notation, we replaced Orb(no) C J^(T>, S*^) with 
^{'D, S"^) itself. Also, we have replaced the isotropy group J^(T>, Doo) with the simpler choice 
J^{V, S0{2)), in analogy with the preceding finite-dimensional treatment. In this case, upon 
fixing no = (0,0, 1), the Lagrange-Poincare reduction with respect to the mechanical connection 
yields the LP Lagrangian fj, n, n, f) on the reduced space 

T Diff™i(D) X {TJ^iV, S^) © T{V)) 3 (r/, ?), n, ri, f) . 

The LP Lagrangian is 

£(?7, ?7, n, fi, f) = i / Wvffi+l / |n|V + ^ / r'^fJ'- F(h o 77"^ V(h o 77-1)/!. 
^ Jv ^ Jv ^ Jv Jv 

At this point, one recalls that the quantity f G J^{T>) in the above Lagrangian is a constant 
parameter, due to the form of the Lagrangian. Thus, it is allowed to set f = for convenience, 
so that the rotational nematic dynamics is regulated by the ordinary Euler-Lagrange equations 
on J^{T>, S"^). Therefore, one obtains the Lagrangian i' : TDiffvoi(^^) x TJ^{V, S*^) — M given by 

f (r7,?7,n,h) : = £(77, ?), n, fi, 0) 

= 11 II^IIV + I / |n|V- / i^(ho7/"\ V(hor/-i))/i 
Jv Jv Jv 

where one can verify directly the invariance of the free energy F under the right action 

(?7, ?7, n, h) ■ (/? = (?7 o y9, ?) o y9, n o y9, fi o ip) 

by simply observing that 

F((h o ip) o [t] o v^)""^, V((h 0^9)0(770 ip)'^y) = F(n o 77^"^, V(n o i]^'^)). 

We now endow the manifold J^{T>, S"^) with the Levi-Civita connection associated to the Rieman- 
nian metric given by integration over V of the natural Riemannian on S*^ C M^. Upon denoting 
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by V/c/t the covariant derivative with respect to the Levi-Civita connection of S^, one writes the 
exphcit form of the Euler-Lagrange equations 



d di' _ _ Q V di' _ V^' _ ^ 
dt di) dr] ' dt dn dn 

on the tangent bundle of Diffvoi(^^) x ^{V, 5*^). The above equations are found for example in 
[15], where they are shown to be equivalent to the Ericksen-Leslie fluid's equations. Of course, 
one can obtain these equations from (5.5), when the initial condition no and are orthogonal. 

5.2.2 Second reduction stage: metamorphosis approach 

The Eulerian form of the fluid equations is found by noting that the Lagrangian i' is invariant 
under the cotangent-lift of the right action {rj, n) i— (77 o y?, n o yj) of the diffeomorphism ip G 
Diffvoi('C)- Thus, one performs the reduction process according to the quotient map 

given by [rj, rj, n, n) H- [u, n, -Djn) := (?) o ry~^, fi o ry~^, n o r/~^) so that the reduced Lagrangian is 

i{u,n,Dtn) = l- I / I An| V " / F{n,Vn)^l. 

^ Jv ^ Jv Jv 

This particular form of Lagrangian reduction is used in [21] to formulate the metamorphosis 
equations in imaging science. See also [12] for a geometric description of metamorphosis reduction. 
Note that we have chosen to denote by Dtn the variable in TnJ^{T>, S^) since the dynamics will 
yield the relation Dtn = {dt + m- V)n. A direct computation shows that the reduced equations 
are given by 

u + u- Vu = —di ( Vn"^ • 7- — I — V», div u 
V 9n^i J 

' YV „\ ^ , , 61 dF 

These equations recover the Ericksen-Leslie fluid equations. 

Remark 5.2 This second reduction step is a standard Lagrangian reduction and not an Euler- 
Poincare reduction as described in Section 2.1, since TJ^{V, S^) is part of the tangent bundle and 
may not be interpreted clS cL parameter space. 

Remark 5.3 (Metamorphosis for parameter-dependent Lagrangians) In order to account 
for compressibility in a natural way, one needs to extend the Lagrangian to depend also on the 
parameter po ^ Den(P), so that 

L(,„,„„) : TDiff (P) X TT{V, 30(3)) ^ M. 

Then, after the first Lagrange- Poincare reduction stage on no = (0, 0, 1), one obtains the reduced 
Lagrangian 

i'p^ : TDiff (P) X {TT{V, S^) © T{V)) -> M. 
32 
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At this point, after setting r = 0, one performs a metamorphosis reduction stage on the above 
Lagrangian according to the invariance property 

i'{r], f], po, n, n) = o 77"^ (po o r]'^) J {rj'^) , n o 77"^ n o 77"^) = i{u, p, n, An) 

where we have used the same notation as in Remark 5.1, although the above reduction process 
differs from that described in 5.1. Therefore, the reduced Lagrangian i is defined such as 

£ : X{V) X Ben{V) x TT{V, S^) R. 

The reduction process just described is an example of a metamorphosis reduction for parameter- 
dependent Lagrangians. Notice that this process can be extended to any type of G-invariant 
Lagrangian 

-.TG xTN ^R. 

where is a manifold, while mo is a parameter belonging to another order parameter space M. 
This method extends the Euler-Poincare approach presented in §2.1. 

5.2.3 Hamilton-Poincare formulation and its Poisson bracket 

Notice that one can Legendre-transform the equations thereby obtaining a Hamiltonian h{m, n, tt), 
on X*oi(X') X T*J^{V,S'^). Indeed, upon introducing the fluid momentum and the conjugate di- 
rector variables 

m = - G XUV) and (n, tt) = (^n, G T*^(P, S') , 

the Hamiltonian functional 

h(?n, n, tt) = (m , n) + (tt , Dtii) — £{u, n, Dtn) 
produces the following Poisson bracket via Legendre transformation 



{f,g}{m,n,Tz) = (m. 



5m'' 6m 



+ ( ./'^ x ,-£jiL(n, 7r))-( £sj^{n,7r)\ (5.10) 

\()(n, tt) sm I \d(n, tt) f>ra I 

where {■, ■}t*j-(ds2) denotes the canonical Poisson bracket on T*T{T>^ S"^), £ denotes Lie deriva- 
tive and the Lie bracket [-, ■] stands for minus the Jacobi-Lie bracket on vector fields, as usual in 
fluid mechanics. The same Hamiltonian structure can also be obtained by a two steps reduction 
from the Hamiltonian : T* ( Diffvoi('C) x J-'{V, 50(3))) — )■ M associated to In this process, 
one proceeds analogously by ignoring the conjugate variable w = 6£/6f (momentum associated 
to f), which is possible because of the special form of the Hamiltonian. 

The Poisson bracket formulation of condensed matter systems (especially liquid crystals) is 
a rather relevant topic in the physics literature; cf e.g. [9, 38]. Here we emphasize that the 
Poisson structures (5.6) and (5.10) arise from two different reductions of the canonical Hamiltonian 
structure on T*(Diff™i(I^) X J^{V,S0{3))) and they both produce the same Ericksen-Leslie 
equations. 
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The Hamiltonian system described by (5.10) is of the general form /i : g* x P — > M, upon 
choosing g = Xvoi(^^) and P = T*J^(V, S'^). However, this construction differs from that treated 
in the first part of this paper and one cannot simply transfer the Lie-Poisson setting discussed 
previously to this infinite-dimensional case. Indeed, the difference resides in the fact that P 
carries its own Poisson structure and thus it is a Poisson manifold by itself. Such a construction 
appears quite often in condensed matter systems, where P = T*J^{V, M) [18], and it also emerges 
in electromagnetic fluid dynamics, where P = T*Q^{V) is the phase space of Maxwell equations 
[17]. This type of Hamiltonian systems arising from an unreduced Hamiltonian on T*G x P has 
been extensively studied in [26], where many interesting properties are presented. For example, 
an interesting consequence of the special type of the bracket (5.10) is that it allows for a Poisson 
isomorphism that eliminates all the terms in the second line of (5.10). More precisely, given 
an equivariant momentum map J : P — )■ g*, there exists a Poisson diffeomorphism {fi,p) ^ 
{fj, + 3{p),p) sending the reduced Poisson structure on (T*G x P)/G = g* x P to the product 
Poisson structure on {T*G/G) x P = g* x P (cf. Proposition 2.2 in [2(:)]). In the case of uniaxial 
nematics, this corresponds to introducing a new variable 

m := m + J(n, tt) , 

so that the new set of variables (m, n, tt) carries the following 'untangled' Poisson bracket: 



{/,^}(m,n,7r) = ( m 



51 6g_ 
5m' Sm 



This bracket produces ordinary Hamilton's equations on T*J^{V, S^) as well as a Lie-Poisson 
equation on Xl^i{V). In the special case of the bracket (5.10) for liquid crystals, one can introduce 
the following momentum map, associated to the cotangent lifted action of G = Diffvoi(I') on 

p = T*T{v,s^y. 

J(n, tt) = Vn^ ■ 77 + V<^ G Xvoi(^^)* (5.12) 

where (n, tt) G T*J-'{V, S"^) and (/? is a variable such that div(J(n, tt)) = 0. Then, the new 
Hamiltonian reads as 

/i (m, n, tt) = ^ I \\m — Vn^ ■ tt — Vip\\'^ ''^ + o~ / l"""!^/^ ~ / -^('^' Vn)yU. 
^ Jv Jv Jv 

This approach has been sometimes referred to as 'untangling' and the Poisson bracket (5.11) 
is called 'untangled Poisson bracket'. See [16, 17, 18] for examples of how the untangling and 
entangling processes are used in the physics of charged fluids and superfiuids. 



5.2.4 The helicity invariant 

Another important property of Hamiltonian systems of the general form h : g* x P ^M. involves 
Casimir functions. Indeed, from the results in [2(:)], it follows immediately that the Poisson bracket 
(5.10) allows for an interesting class of Casimir functions. Indeed, given an equivariant momentum 
map J : P — 7- 0* and any Casimir function C(/i) for the Lie-Poisson bracket on g*, the function 

C(/i,p) = C(/i + J(p)) (5.13) 
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is a Casimir for the reduced Poisson bracket on {T*G x P)/G = q* x P (cf. Corollary 2.3 in 
[26]). In the case of liquid crystals, the momentum map (5.12) can be used to produce an explicit 
expression for the helicity of nematic liquid crystals. Indeed, since it is well known that the 
helicity Jif{m) = (curl m,m) is a Casimir for the Lie-Poisson bracket on Xvoi(^^)*, then direct 
substitution of the momentum map J(n, tt) in formula (5.13) yields a Casimir for the Poisson 
bracket (5.10). This Casimir is explicitly written as the following helicity functional for nematic 
liquid crystals 

^(m, n, 7r)= / (m + Vn^ ■ tt) ■ curl (m + Vn^ ■ tt) . (5.14) 
Jv 

This helicity invariant can also be written as a Casimir for the Lie-Poisson bracket (5.6). Indeed, 
one can pull back using the mapping (m, ^t, n) H- (m, n, /x x n), which is suggested by the 
isomorphism (4.18) (with w = 0) holding for the dynamics of the single nematic molecule. The 
resulting expression for the helicity is then 

Jif{m, /X, n) = / (m + /X ■ (n X Vn)) ■ curl (m + /x ■ (n x Vn)) . (5.15) 
Jv 

More rigorously, the Casimir property can be verified by restricting the Lie-Poisson construction 
in (5.6) to the submanifold = {(m, /x,n) | /j, ■ n = 0}, which is analogous to setting w = 
in (4.18) for a single nematic molecule. Then, the following arguments show that ^ can be 
endowed with the Poisson structure (5.6), thereby making it into a Poisson submanifold of W = 
^voi(^^)* X J^{T>,so{3))* X T{V,S'^). Let {m{t) , ^j,(t) , n{t)) be a solution of Hamilton's equation 
/ = {/, h} on W relative to an arbitrary Hamiltonian h. Since dt{fi ■ n) + (n • V)(/x • n) = 0, 
then we obtain {n{t) ■ n(t)) = (^^(0) ■ n(0)) o r]^-^, where we have denoted by r]t the flow of u. 
Consequently, any Hamiltonian vector field on W restricted to is tangent to and this shows 
that ^ is a quasi Poisson submanifold. Thus, there is a unique Poisson structure on ^ making 
it into a Poisson submanifold, see Prop 4.1.23 in [3G]. It is readily seen that this Poisson structure 
has the same expression (5.6). 

Theorem 5.4 The expression 

M'[m, /X, n) = / (m + ■ (n X Vn)) • curl (m + ■ (n x Vn)) . (5.16) 
Jv 

is a Casimir function on the Poisson manifold 0^ = {(m,/x, n) | /x ■ n = 0}, endowed with the 
bracket (5.6). 

Proof. We shall show that for a solution (m(t), ^i(t), n(t)) of an arbitrary Hamiltonian system 
on we have 
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where we have defined the differential one form C := m + fi ■ {n x Vn). A direct calculation shows 
that 

9 \^ 5h ^ 6h , fSh Sh 

— + £sh C = -/x- — + Vn-- gradp + ^ x M + ^ 

ot j:;^ J on on \o/^ on 

J J \ 

— /X • h Vn ■ gradp — • v — ■ n n — V — 



Sh 


+ Vn- 


Sh 




U ±1. 


(( 


Sh 

— X n 
Sfi 


)' 


5h 


+ Vn- 


Sh 


6fi 


Sn 


6h 


+ Vn- 


Sh 


6fi 


Sn 



Sfj, J Sfj, J 

= — dp 

where we have used the Jacobi identity for double cross products, and the properties |np = 1 and 
/X - n = n - Vn = 0, respectively. Thus, upon writing the helicity as = C A dC, one obtains 

dt {C A dC) = -dp AdC- div 1 (C A dC) 1^1 d=^a; 

\ om / 

so that = J^C A dC = const along any Hamiltonian flow on J^. M 
We notice that a simple consequence of the above theorem is the 

Corollary 5.5 With the initial condition /Xq - no = 0, the circulation theorem for the equations 
of liquid crystals reads as 

— (b (m + /X - (n X Vn)) = , 

which follows directly from the property {pt + £sh/5m) C = —dp. Similar arguments on the 
relation between the Kelvin-Noether circulation and the helicity invariant also hold for superfluid 
dynamics [20]. 

Remark 5.6 (Two dimensional flows) It is important to notice that formula (5.13) provides 
a whole class of Casimir functions for 2D flows. Indeed, it is well known that, in two dimensions, 
any function $(0;) of the vorticity u = z ■ curl(m) yields a Casimir C = J $(a;) and thus, by 
(5.13), any function $(a; + z - curl(Vn - tt)) produces a new Casimir 

C = [ $(c^ + {7r„n"})= / ^{u + eabc{cT'n',n''}) 
Jv Jv 

for 2D flows of uniaxial nematics. Here, the operation {-, -} is the canonical Poisson bracket in 
the planar coordinates (x, y). This observation may have important consequences in terms of the 
stability properties of liquid crystals, following the approach in [23]. 



36 



5.3 Reduction schemes for uniaxial nematic liquid crystals 

Starting from the same Lagrangian L^^ on the tangent bundle of Diffvoi('C') x J'iV, S0{3)) we 
have obtained the Ericksen-Leshe fluid equations in two different ways. For completeness, this 
section considers nonvanishing r. 

First, by considering as symmetry group the semidirect product Diffvoi(I^) © 50(3)), 
we obtain the equations of motion by Euler-Poincare reduction, following the general approach 
described in §2.1. The reduced space is given by (Xvoi('C) (S) so(3))) x J^(V,S^) and the 
reduction map reads 

iv, V, X, X, no) I — y {{r], ?), X, x){v, x)"\ iv, x)^o) 

= (?) o 77-1, (xx~^) o ?7~\ (xno) o 7]-^) 
=■■ {u,i^,n) 

according to the general formula (2.1). This approach is consistent with that given in [13]. 

Second, we use Lagrangian reduction by stages (see [^)]), and obtain the motion equations on 
the reduced space Xvoi('C) x TJ^{V, S^) x J^(T>). First we apply Lagrange-Poincare reduction for 
the internal structure, with respect to the mechanical connection (here we choose no = (0,0, 1)). 
We get 

(??, ?7, X, X, no) I — > (?7,?7,xno, {xx'^) x ix^o) , ixi^o) ■ XX'^) 
=: (?7, ri,n,u X n, u n) =: (77, fj, n, n, f). 

Then, we use Lagrangian reduction with respect to the fluid variables and have 

(?7, ?7, n, n, f) 1 — > (rj o 77"^, n o 77^^, n o ri~^,f o r]"^) =: (m, n, DtU, r) . 

This approach is consistent with that given in [15]. Note the following important relation between 
the intermediate and fully reduced variables. We have 

n = xno, n = (^hq) o r^^^ = n o rj'^ 

f = (xno) -xx"^ = i>-n, r = u-n = f o r]~^ . 

One can pass from one approach to the other by using the same transformations (4.13) and 
(4.14) as for the single particle case, except that now they are applied to vector valued functions 
instead of vectors. This clarifies the link between the approaches given in [15] and [13]. The 
diagram below clarifies the various reduction processes involved in the hydrodynamics of nematic 
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liquid crystals. 



TDiff™i(^^) X TJ-(P, ^0(3)) 



TDiff™i(r') X {TJ^{V,S'')®J'{V)) 



(X™i(P)(s)^(I?,so(3))) X Xvoi(2^) X (TJ-(I),52) ©^(P)) 

(r) o ?7~\ i> o ?7~\ n o 77"^) =: [u, u, n) (?) o 77"^ n o 77"^ fi o 77"^ f o 77"^) =: (n, n, An, r) 

6 Conclusions 

This paper has developed symmetry-reduction methods for systems with broken symmetry, whose 
main example consists of nematic particles in liquid crystals (ordinary or biaxial). After extending 
the Euler-Poincare theory, this has been applied to the case of a transitive group action on an 
order parameter manifold, thereby showing how the latter carries its natural coset structure 
typically appearing in symmetry breaking. This setting has been applied to nematic particles 
(5*0 (3) /O (2)), biaxial nematics {SO{3)/D2) and V^-shaped molecules (5*0 (3) /Z2), and its general 
validity can be used to reproduce dynamics on any order parameter space. Moreover, the Euler- 
Poincare and Lie-Poisson dynamics for nematic and biaxial systems have been formulated also in 
terms of the alignment tensors, which are widely used in the Landau-de Gennes theory. 

As a further step, the Lagrangian reduction technique (Lagrange-Poincare reduction) has 
been applied directly to the space (TG)/Go, where Go C G is the isotropy group determining 
the "breaking symmetry" , which is left in the system when the full G-symmetry is broken. This 
technique has showed how the reduced configuration space may be identified with the coset G/Gq 
(order parameter space), provided the adjoint bundle term in can be set to zero. This is 
indeed the case for nematic particles, so that the resulting dynamics is given by the Ericksen- 
Leslie equation on S'0(3)/0(2) for a single director. The choice of the mechanical connection 
enabled us to explain why Ericksen-Leslie dynamics can be considered as ordinary Hamiltonian 
dynamics on the projective plane, upon setting n - u = for consistency with the rod-like nature 
of nematic particles. Moreover, the mechanical connection produced explicit transformations that 
relate directly the Euler-Poincare and the Lagrange-Poincare approaches. The Lagrange-Poincare 
approach restricts to the case when the isotropy subgroup Gq C G is a Lie group with dimGo > 1, 
which excludes the case when Gq is discrete. It is then an interesting open question whether it 
is possible to extend the Lagrange-Poincare reduction process (TG) / Gq to account for a discrete 
symmetry group Gq. 

All the above techniques have been extended at the continuum level, to produce also hydro- 
dynamic models, thereby recovering the well known Ericksen-Leslie equation for the dynamics 
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of the director field. As a result of this extension, all the geometric features of the microscopic 
particle dynamics are naturally transfered at the macroscopic fluid level without substantial mod- 
ifications. Interestingly enough, the application of Lagrangian reduction to the hydrodynamics 
of nematic liquid crystals produces a "generalized Lie-Poisson" Hamiltonian structure on q* x P, 
involving a Lie algebra q and a Poisson manifold P. This type of Hamiltonian structure was 
discovered in [26] and it appears very frequently in hydrodynamic physical models as well as 
in other contexts, such as imaging sciences. The last part of this paper applied the theory on 
this generalized Lie-Poisson construction to produce new Poisson bracket structures, as well as 
new Casimir functions for liquid crystal dynamics. In particular, two explicit expressions for the 
helicity of uniaxial nematic liquid crystals have been presented in (5.14) and (5.15). 
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